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G^ Abstract 
O 

Recent Nuclear Magnetic Resonance and elastic neutron scattering experiments have revealed conclusively the presence of static 
incommensurate magnetism in the field-induced B phase of CeCoIns. We analyze the NMR data assuming the hyperfine coupling 
to the In(2) nuclei is anisotropic and simulate the spectra for several different magnetic structures. The NMR data are consistent 



with ordered Ce moments along the [001] direction, but are relatively insensitive to the direction of the incommensurate wavevector. 
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The heavy fermion superconductor CeCoIng exhibits a rich 
spectrum of strongly correlated electron behavior. This uncon- 
ventional d-wave superconductor exhibits non-Fermi liquid be- 
havior associated with proximity to a proposed quantum critical 
point, as well as a new thermodynamic phase (B phase) that ex- 
ists only within the superconducting phase near i/,2 fuH- Ini- 
tially this B phase was identified as the elusive Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) superconducting phase first pre- 
dicted to exist in Pauli-limited superconductors over 40 years 
iigo m 0. 13. Si- However, recent NMR work identified the 
presence of incommensurate antiferromagnetic order in the B 
phase in contrast to the standard predictions for the FFLO phase 
OlslISI]. Signatures of magnetism were also seen in other NMR 
experiments ifioi [Till . Despite initial arguments to the con- 
trary II12II . recent neutron scattering results by Kenzelmann and 
Coworkers now provide conclusive proof for long-range static 
incommensurate antiferromagnetic order lll3ll . 

The original NMR work measured the spectrum of the In(2) 
sites and proposed a candidate magnetic structure in which both 
the ordered Ce moments and the incommensurate wavevector 
are parallel to the applied field (along [100]). However, the 
neutron scattering experiments found that when the field was 
applied along [110] the moments lie along [001] and the in- 
commensuration along [110]. Recently we have re-analyzed the 
NMR data and shown that by allowing for anisotropic hyperfine 
coupling between the Ce spins and the ' '^In nuclei, rather than 
simple isotropic coupling as originally assumed in [7], then the 
NMR spectra are fully consistent with moments along [001] 
lOI- Here, we show that by using hyperfine parameters con- 
sistent with experimental measurements of the Knight shift, we 
find quantitative agreement for the "''In NMR spectra and put 
constraints on the incommensuration wavevector. 




Figure 1 : (Color online) The unit cell of CeCoInj . The Ce atoms (yellow) sit 
at the eight corners. The In(l) atoms sit in the center of the top and bottom 
faces (orange). The Co atoms are grey and the In(2) atoms are green. For the 
field oriented in the ab plane, there are two inequivalent In(2) atoms, depending 
on whether the field is parallel, In(2a), or perpendicular, In(2b), to the unit cell 
face. 



2. Incommensurate Antiferromagnetism 

The incommensurate magnetic structure is given by S(r) = 
Socos[(Qo + Qi) ■ r], where So is the ordered moment, with 
commensurate, Qo = {7r/a,7r/a,n/c), and incommensurate, 
Q, - 7T/a{6x,6y,Q), antiferromagnetic wavevectors. Kenzel- 
mann et al. found 5^ - Sy - 0.12 and - O.lSyU^ along 
[001]. The ordered moments create a static field at the "^In and 
^^Co nuclear sites through the hyperfine coupling. This field 
cancels by symmetry at the In(l) and Co sites, but the NMR 
spectra show a large hyperfine field at the In(2a) site that is 
aligned either parallel or antiparallel to the applied field (along 
[100]). The hyperfine coupling to the In(2a) site is given by 
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Figure 2: (Color online) The static hyperfine field at the In(2a) site along the 
[100] (blue), [010] (purple) and [001] (brown) directions as a function of lat- 
tice position along [100] for So = 0. 15/jb and 6y = Sy = 0.121, with lattice 
parameter a, respectively. The solid lines show the modulation of the staggered 
magnetization (Eq.|2). 



"H = 'Eienn I ' B ' S(R,), where the sum is over the two nearest 
neighbor Ce spins (along [100]). The coupling tensor is given 

by B = Bi,oI + Bdip, 



(1) 
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The constants Bis,,, B^ip and 6- can be determined by fits to the 
Knight shift and are given by: Bis,, = ^rf/p = 7.1 kOe/fis and 
6. - 29° jT3]. In the B phase, the static hyperfine field at the 
In(2a) site varies spatially with 

Hhfila) = SoBdip[3 sin(26i,)cos(-^)a 



((l-cos(20,))sin(^))c]. (2) 



In order to make detailed comparisons with experimental 
spectra, we have calculated the hyperfine field and the reso- 
nance frequency, / = y|Ho + H/,/|, where yHo || [100] ~ 118.3 
MHz, for several different values of 6x and for a lattice of 100 
X 100 unit cells. Representative cuts in real space are shown in 
Fig. |2] Spectra were determined by calculating the histogram 
of frequencies, shown in Fig. |3] As seen in Fig. |2] the hyper- 
fine field at the In(2a) site has components either parallel or an- 
tiparallel to the applied field leading to the broad double-peak 
spectra seen in Fig. [3] However, the size of the ordered mo- 
ment measured by neutrons. So = 0.1 Syu^ implies a hyperfine 
field larger than observed by NMR. The NMR data are more 
consistent with an ordered moment of O.QlfiB- Using Eq. |2]and 
the measured value of the modulus of the hyperfine field of 1.3 
kOe, we find that the ordered moment and incommensuration 
are related by 6^ = 5v0 cos"'(5q/5o), where Sm = 0.64 and 
= 0.07 fiB (shown in Fig. Elf)). 

As seen in Fig. |3la) -|3le), the measured spectra is relatively 
insensitive to the direction of the incommensurate wavevector 



Q, . Subtle changes in the calculated spectra arise because of the 
component of H/,y (2fl) along [001] gives rise to small frequency 
shifts. 

The difference between the magnitude of the ordered mo- 
ment measured by neutrons versus that measured by NMR may 
be due to either (i) changes in the magnetic structure as a func- 
tion of the orientation of the field, (ii) errors associated with 
the measured values of the hyperfine coupling, (iii) dynamic 
fluctuations of the magnetic structure, or (iv) neglect of cou- 
pling to conduction electrons. There are two types of hyperfine 
couplings in heavy-fermion compounds, one to the conduction 
electron spins, S,-, and another to the local moments, S/ ITsI] . 
Here, we ignored the coupling to Sc, since the static hyperfine 
field presumably arises from the ordered Ce moments in the 
B phase. Including a hyperfine coupling to the Sc spins may 
account for the difference in the ordered moments and the in- 
commensuration indicative of a Fermi surface instability. 

In the original measurements of the In(2a) NMR, the broad 
spectrum only emerged below ~ 100 mK (at 11.1 T where 
Tc = 470 mK and To = 290 mK, where To is the transition 
temperature to the B phase). Above ~ 100 mK the In(2a) spec- 
trum disappeared. The reason for this behavior is likely due to 
dynamic effects, in which the spin-spin relaxation rate ' is 
excessively fast due to the dynamics of the magnetic structure. 
For sufficiently low temperatures where the magnetic structure 
becomes static and ' is small enough the broad spectrum 
emerges. However, when the time scale of the fluctuations is 
of the inverse of the linewidth (~ 2.8 kHz in this case) the 
splitting of the NMR spectrum can be less than the full static 
linewidth. It is possible that for temperatures T «; 60 mK, the 
NMR linewidth would broaden to reflect a larger ordered mo- 
ment. 

3. Superconductivity 

An outstanding question regarding the nature of the B phase 
is the possible modulation of the superconducting order param- 
eter A. Using Ginzburg-Landau theory, Kenzelmann et al. pro- 
posed a specific coupling between A, the antiferromagnetic or- 
der parameter Mq, and the magnetic field H to account for the 
phase diagram [13]. They postulated that at large fields an in- 
duced subdominant order parameter nucleates with finite mo- 
mentum, Aq, in addition to the dominant, uniform d-wave order 
parameter Aq. In real space, this scenario suggests that the su- 
perconductivity and antiferromagnetism are out of phase, such 
that the extrema of magnetization Mq are located at the nodes 
of Aq. Consequently, Aq should vary at atomic length scales, 
which is physically unreasonable considering the stiffness of 
the condensate, 27r/|Q| w 2a «; ^q. 

The NMR spectra of the In(l) site, where the hyperfine field 
cancels, show no broadening that would be associated with a 
spatial variation of the Knight shift K due to real-space nodes 
iLZlj. is proportional to the spin susceptibility, which vanishes 
for T «: A/ks [16]. If A varied spatially and vanished at real- 
space nodes, then K would vary as well leading to a variation 
of the resonance frequency and a broad resonance |@|. An alter- 
native scenario is that A retains a large uniform component, but 
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Figure 3: (Color online) Simulated spectra and experimental data at 11.1 T and 50 mK (solid points, reproduced from 01) for (a) {S_^,Sy) = (0,0.121) and 
So = O.OVps, (b) (6^,6y) = (0.121,0.121) and So = 0.07/Js, (c) {S^,6y) = (0.121,0.121) and So = O.lSpjj, (d) (Sj^Sy) = (0.498,0.498) and So = O.lO^jg, (e) 
i^jc, ^y) = (0.681, 0.681) and So = 0. 15yUg. (f) The locus of points (So, 6x) that are consistent with experiment. 
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has an induced component with modulation Q,: A = Aq + Aq,, 
with 2;r/|Q,| > ^o- Our data suggest that Aq is large and discon- 
tinuous across the first order normal-to-superconducting transi- 
tion at Tc ~ 470 mK at 1 1 . 1 T. The narrow linewidth of the In(l) 
(~ 33 kHz at 100 mK) suggests that \dK\/(K(T,)-K(0)) < 0.07, 
where 6K is the spatial modulation of the Knight shift due to 
Aq_ . The absence of any Knight shift modulation at the In(l) 
preclude any real-space nodes in the order parameter A and are 
at odds with an exchange-spUt Fermi surface resulting in an 
FFLO state. 



4. Discussion 

The NMR data suggest that the B phase of CeCoIns con- 
sists of incommensurate antiferromagnetism coexisting with d- 
wave superconductivity in the vortex state. This material ex- 
hibits a large Fermi surface consistent with fully hybridized f 
electrons [17], thus the presence of ordered local moments with 
magnitude 0.07 fig in the superconducting state remains unclear 
Possible explanations include modulated hybridization lUsIl or a 
two-component scenario [19]. The emergence of intrinsic spa- 
tial inhomogeneity in pristine undoped materials under these 
conditions represents a fascinating consequence of proximity 
to a first-order quantum critical phase transition, which may or 
may not be universal [20, 21]. 

We thank R. Movshovich, V. Mitrovic, A. Balatsky, L. 
Boulaevskii and M. Kenzelmann for valuable discussions and 
sharing their results. Work at Los Alamos National Laboratory 
was performed under the auspices of the US Department of En- 
ergy under grant no. DE-AC52-06NA25396. 
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